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Abstract. We utilize a 3D hybrid hydro+micro model for a comprehensive and 
consistent description of soft and hard particle production in ultra-relativistic heavy- 
ion collisions at RHIC. In the soft sector we focus on the dynamics of (multi-)strange 
baryons, where a clear strangeness dependence of their collision rates and freeze-out is 
observed. In the hard sector we study the radiative energy loss of hard partons in a soft 
medium in the multiple soft scattering approximation. While the nuclear suppression 
factor Raa does not reflect the high quality of the medium description (except in a 
reduced systematic uncertainty in extracting the quenching power of the medium), the 
hydrodynamical model also allows to study different centralities and in particular the 
angular variation of Raa with respect to the reaction plane, allowing for a controlled 
variation of the in-medium path-length. 



Experiments at the Relativistic Heavy Ion Collider (RHIC) have established a 
significant suppression of high-p^ hadrons produced in central A-l-A collisions compared 
to those produced in peripheral A-f A or binary scaled p-|-p reactions, indicating a 
strong nuclear medium effect [H |2], commonly referred to as jet- quenching. Within 
the framework of perturbative QCD, the leading process of energy loss of a fast parton 
is gluon radiation induced by multiple soft collisions of the leading parton or the radiated 
gluon with color charges in the quasi-thermal medium [3], HJ [5] . 

Over the past two years, a large amount of jet-quenching related experimental data 
has become available including but not limited to the nuclear modification factor Raa, 
the elliptic flow V2 at high px (as a measure of the azimuthal anisotropy of the jet 
cross section) and a whole array of high hadron-hadron correlations. Computations 
of such jet modifications have acquired a certain level of sophistication regarding the 
incorporation of the partonic processes involved. However, most of these calculations 
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have been utilizing over-simplified models for the underlying soft medium, e.g. assuming 
a simple density distribution and its variation with time. Even in more elaborate 
setups, most jet quenching calculations assume merely a one- or two-dimensional Bjorken 
expansion. 

The availability of a three-dimensional hydrodynamic evolution code |6j and related 
hybrid approaches allow for a much more detailed study of jet interactions in a 
longitudinally and transversely expanding medium. The variation of the gluon density 
in these approaches is very different from that in a simple Bjorken expansion. A previous 
calculation in this direction [TJ [8] estimated the effects of 3-D expansion on the Raa- 
However, this approach treated the energy loss of jets in a rather simplified and heuristic 
manner. Here, we shall perform a detailed investigation of the modification of jets 
in a three dimensionally expanding medium within the BDMPS formalism utilizing 
quenching weights as described in [9|. In addition, 3-D hydrodynamic and hybrid models 
have been very successful in describing the majority of features of soft particle production 
at RHIC (with HBT interferometry being the sole exception) - this we shall utilize in 
order to determine the medium properties for the jet-quenching calculation. 

In this write-up, we utilize a state-of-the-art fully 3-D hybrid hydro-|-micro 
transport model [6]. The model employs relativistic 3D-hydrodynamics for the early, 
dense, deconfined stage of the reaction and a microscopic non-equilibrium model for 
the later hadronic stage where the equilibrium assumptions are not valid anymore. It 
is capable of self-consistently calculating the freezeout of the hadronic system, while 
accounting for the collective flow on the hadronization hypersurface generated by the 
QGP expansion. The initial conditions of the hydrodynamic calculation are tuned to 
describe the hadronic data in the soft sector, such as hadron yields, spectra, rapidity- 
distributions as well as radial and elliptic flow. Note that while the modeling of the 
hadronic stage is of paramount importance for the proper description of the medium in 
the soft sector, it does not contribute to the jet energy loss. 

In the left frame of Fig. [1] we analyze the Pt spectra of multistrange baryons. 
Our results show good agreement with experimental data for A, S, Q from the STAR 
collaboration [TOl |Tl]. Recent experimental results suggest that at thermal freezeout 
multistrange baryons exhibit less transverse flow and a higher temperature closer to 
the chemical freezeout temperature compared to non- or single-strange baryons [TOl [H] . 
This behavior can be understood in terms of the flavor dependence of the hadronic 
cross section, which decreases with increasing strangeness content of the hadron. The 
reduced cross section of multi-strange baryons leads to a decoupling from the hadronic 
medium at an earlier stage of the reaction, allowing them to provide information on the 
properties of the hadronizing QGP less distorted by hadronic final state interactions 
[T2| [T3l [H]. In microscopic calculations the early decoupling will manifest itself via a 
reduced number of collisions for the respective hadron species. It should be noted that 
the analogous behavior has already been observed in experiments at the CERN-SPS 

[El [161 [mini [19]. 

The pseudo-rapidity dependence of the number of hadronic rescatterings for 
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Figure 1. Left: Pt spectra of (niulti-)strange baryons calculated in a hybrid 
hydro+micro approach compared to STAR data. Right: pseudo-rapidity distribution 
of the number of rescatterings for baryons of varying strange quark content. 



different baryon species can be used to corroborate the above findings and is analyzed 
in the right frame of Fig. [H which shows the number of colhsions of p, A, S and as a 
function of ?7 at 6 = 2.4, 4.5 and 6.3 fm. The distributions appear to be similar to that of 
the particle yield psuedorapidity distribution. At midrapidity we find a plateau region 
extending from rj = —3 to 3, followed by a steep drop-off to forward and backward 
rapidities. The flavor dependence of the average collision numbers is again clearly seen, 
even though we would hke to point out that the shapes of the different distributions are 
very similar. The large plateau region indicates the rapidity domain in which interacting 
matter can be found and in which the application of thermodynamic concepts is viable. 

Having determined the properties and dynamics of the soft sector, we can now 
utilize the time-evolution of the medium provided by our model for the calculation of 
jet energy-loss. Thus, our calculation significantly reduces the systematic uncertainties 
usually associated with the medium parametrization and allows for a precision 
calculation of all effects associated with hard probe - medium interactions. Our 
calculation follows the BDMPS formalism for radiative energy loss [20] using quenching 
weights as introduced by Salgado and Wiedemann [211 IS] • 

In [22l l23j it has been shown that Raa for central collisions only constrains a scale, 
but not the detailed functional form of energy loss probability distribution {P{AE))tj^^- 
In the approach outlined above, this is manifest in the parameter K in the expression 
for the local transport coefficient q{C,)'- 

m =K-2. (1) 

which then was adjusted to the data in central collisions. We illustrate in the left 
frame of Fig. [2] that three different dynamical models, a 2D hydrodynamical evolution 
[24] , the 3D hydrodynamical evolution outlined above [6] and a parametrized fireball 
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Figure 2. Left: Raa for central collisions as calculated in three different models 
for the medium evolution with the overall quenching power scale K adjusted to data. 
Right: pt dependence of Raa in plane (solid) and out of plane (dashed) emission at 
different values of impact parameter b. 



evolution [25] give almost equal descriptions of Raa once the scale parameter is adjusted, 
albeit they require different values of K (the chief reason for this being the different 
longitudinal dynamics). The ±50% spread in the values of if for the different models of 
the medium is a measure for the systematic error inherent in the tomographic analysis 
of jet energy-loss via the nuclear modification function Raa- 

However, one may gain predictive power in going to collisions at finite impact 
parameter b. In particular, the dependence of Raa for non-central collisions 
constitutes a systematic variation of path-length within a system with fixed overall 
scale. The average path-length is expected to be smaller for a parton emitted in plane 
as compared to one emitted out of plane, and hence Raa is expected to be larger at 
= than at = 7r/2 with the difference in Raa between these angles increasing with 
the initial asymmetry (and hence b). Using a simple model for the time-evolution of 
the medium and collective flow effects, it has been shown in [26] that the dependence 
of Raa is quite sensitive to the initial gluon density distribution and temporal evolution 
of the medium. 

Utilizing the previously discussed 3-D RFD model [6] , we show the pt dependence 
of Raa for emission in plane and out of plane at three different impact parameters b 
in the right frame of Fig. [2j As expected, Raa grows for more peripheral collisions as 
there is less soft matter produced to induce energy loss. Moreover, there is a smooth 
angular variation of Raa observed, reflecting the underlying medium asymmetry. The 
difference between in-plane and out of plane emission grows with impact parameter, at 
b = 2.4 fm there is hardly angular variation whereas at 7.5 fm differences are of order 
20%. 

In summary, we have utilized a 3D hybrid hydro+micro model for a comprehensive 
and consistent description of soft and hard particle production in ultra-relativistic heavy- 
ion collisions at RHIC. In the soft sector we have focused on the dynamics of (multi- 
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) strange baryons, where a clear strangeness dependence of their colhsion rates and 
freeze-out is observed. In the hard sector we have studied the radiative energy loss 
of hard partons in a soft medium in the multiple soft scattering approximation. Our 
analysis should be seen as the starting point for a comprehensive study of probe-medium 
interactions treating the hard and soft sector on equal footing. 
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